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Анализируются методы расчета температурных и потенциальных полей, а также вольт-амперных характеристик многоэлемент- 
ных термоэмиссионных электрогенерирующих сборок в вакуумном, диффузионном и разрядном режимах работы. Приведено 
описание алгоритма решения самосогласованной двумерной задачи тепло- и электропроводности со сложными граничными 
условиями, позволяющее получить зависимости генерируемой мощности от удельного тепловыделения, давления пара цезия, 
сопротивления нагрузки, входной температуры и расхода теплоносителя. 


Для расчета вольт-амперных характеристик 
(ВАХ) многоэлементных электрогенерирующих 
сборок (ЭГС) существует несколько подходов [1—7]. 

Первый из них [1] состоит в расчете серий тем- 
пературных и потенциальных полей и интеграль- 
ных ВАХ отдельных электрогенерирующих эл- 
ементов (ЭГЭ). Учитывая, что по элементам при их 
последовательном соединении проходит одинако- 
вый ток, далее используется метод суммирования 
напряжений отдельных ЭГЭ для задаваемых вели- 
чин тока. Недостатком такого метода является 
необходимость расчета большого количества ВАХ 
отдельных элементов, особенно при неравномер- 
ном тепловыделении и применении геометриче- 
ского профилирования ЭГС. 

Второй способ основан на совместном решении 
уравнений теплопроводности и электропроводно- 
сти сборки в целом [2-7]. К достоинствам второго 
подхода относятся автоматический учет неравно- 
мерного распределения тепловыделения, гранич- 
ных условий между элементами, возможность мо- 
делирования возникновения и развития аварийных 
ситуаций [3], т.к. ЭГС рассматривается как единое 
целое. 

При выводе термоэмиссионного реактора-пре- 
образователя на проектную мощность ЭГС работа- 
ют в различных режимах: вакуумный (квазиваку- 
умный) и цезиевый (плазменные). Поэтому для со- 
провождения испытаний ЭГС разрабатываемые 
методы расчета ее характеристик должны учиты- 
вать следующие особенности [8]. 

Во-первых, на этапе обезгаживания эмиттер- 
ных оболочек при давлениях в полостях межэлек- 
тродного зазора (МЭЗ) порядка 10 4 ... 10 5 Па и тем- 
пературах эмиттера примерно 1800... 1900 К, а то- 
плива 1900. ..2000 К, производится съем вакуумных 
ВАХ ЭГС. По характерным точкам и наклону ваку- 
умных ВАХ определяются неизмеряемые параме- 
тры, а именно вакуумные работы выхода эмиттера 
и коллектора, температура эмиттера, сопротивле- 
ния утечек тока через дистанционаторы и коллек- 
торную изоляцию, а также контролируются ваку- 
умные условия в МЭЗ. 

Во-вторых, после окончания термовакуумной 
подготовки производится перевод ЭГС в цезиевый 


режим работы с постепенным увеличением давле- 
ния паров цезия в МЭЗ. При малых значениях да- 
вления паров цезия, соответствующих температу- 
рам примерно до 500 К, ЭГС работает в прямопро- 
летном либо диффузионном режиме. В этих режи- 
мах также производится съем ВАХ с целью оценки 
значений внутренних параметров ЭГС. 

В-третьих, при повышении температуры пара 
цезия примерно до 520 К в МЭЗ отдельных ЭГЭ во 
время съема ВАХ происходят микроподжиги раз- 
ряда, поэтому на ВАХ наблюдаются так называ- 
емые ступеньки - скачки тока и напряжения. По 
количеству ступенек оценивают число ЭГЭ, пере- 
шедших в дуговой, основной режим работы. 

Наконец, в связи с отсутствием возможности 
экспериментального измерения внутренних пара- 
метров (работы выхода и температуры электродов, 
межэлектродный зазор и т.д.), определяющих энер- 
гетические характеристики ЭГС, одним из эффек- 
тивных методов их оценки в настоящее время счи- 
тается метод сравнения экспериментальных ВАХ с 
расчетными. Причем, по расчетным характеристи- 
кам определяются не только усредненные параме- 
тры, но и распределения температур, потенциалов 
и токов, которые позволяют выявить наиболее тер- 
мо- и электронапряженные участки конструкции. 

Поэтому проблема разработки эффективных 
алгоритмов расчета характеристик ЭГС, позволяю- 
щих имитировать экспериментальные исследова- 
ния в различных режимах работы (вакуумный, 
диффузионный, переходный, дуговой, переключе- 
ние нагрузки, съем ВАХ с различной скоростью и 
т.д.) является актуальной. Повышение эффектив- 
ности алгоритмов может быть достигнуто за счет 
использования особенностей работы ЭГС в раз- 
личных режимах. 

Например, существенной особенностью ваку- 
умного и диффузионного режимов работы являет- 
ся низкий уровень плотности генерируемого тока, 
примерно 0... 10 6 и 0.. .0,2 А/см 2 , соответственно. 
Поэтому электронная составляющая теплового по- 
тока с эмиттера практически не оказывает влияния 
на температурное поле ЭГС. Кроме того, выделе- 
ние тепла при прохождении тока по электродам и 
коммутационным перемычкам незначительно, а 
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электроды можно считать изопотенциальными. То 
есть температурное поле электродов во время съе- 
ма ВАХ не изменяется. Эти обстоятельства позво- 
ляют раздельно анализировать тепловые и элек- 
трические процессы в ЭГС и применить для расче- 
та вакуумных и диффузионных ВАХ следующий 
подход. По заданным внешним параметрам (уро- 
вень и профиль тепловыделения, условия охлажде- 
ния, давление пара цезия) рассчитывается темпе- 
ратурное поле (для диффузионного режима с уче- 
том теплопроводности паров цезия). По получен- 
ному температурному полю электродов и заданной 
величине тока находятся напряжения отдельных 
элементов с учетом возможных утечек тока между 
электродами. Полное напряжение ЭГС определя- 
ется методом суммирования напряжений отдель- 
ных ЭГЭ. 

В диффузионном режиме с более высокими 
плотностями генерируемого тока, а особенно в ду- 
говом, существенную роль в тепловом балансе 
электродов начинает играть электронное охлажде- 
ние и джоул ево тепловыделение, а также проявля- 
ется неизопотенциальность электродов. Проблема 
расчета ВАХ ЭГС, соответствующих переходу от 
диффузионного режима к дуговому, рассматрива- 
ется в работе [4] . В ней показана принципиальная 
возможность моделирования ступенчатых характе- 
ристик, связанных с поджигами разряда в отдель- 
ных элементах, причем результаты моделирования 
подтверждаются экспериментальными данными. 

Таким образом, методика расчета ЭГС в ваку- 
умном и диффузионном режимах сводится к реше- 
нию уравнений теплопроводности с целью получе- 
ния температурных полей электродов. Далее на ос- 
нове использования вакуумных либо диффузион- 
ных ВАХ термоэмиссионного преобразователя 
(ТЭП) производится расчет ВАХ ЭГС. Для дугово- 
го режима при расчетах статических ВАХ и пере- 
ходных процессов необходимо учитывать взаимное 
влияние тепловых и электрических процессов. Ис- 
ключением является расчет изотермических ВАХ с 
высокими скоростями, когда за время их съема 
температурное поле не успевает измениться. 

Для расчета статических и динамических харак- 
теристик ЭГС необходимо найти самосогласован- 
ное решение уравнений, приведенных в [9] с соот- 
ветствующими граничными и начальными усло- 
виями при внесении различных возмущений в ма- 
тематическую модель сборки. 

Представление математической модели ЭГС в 
виде нестационарной тепловой и стационарной 
электрической позволяет разбить сложную задачу на 
две более простых. Первая задача сводится к нахож- 
дению распределения потенциалов по электродам 
сборки при известном распределении их температур 
и заданному сопротивлению нагрузки. Вторая задача 
заключается в определении решения уравнений те- 
плопроводности с соответствующими граничными 
условиями с учетом полученного распределения то- 
ков и потенциалов по высоте ЭГС. При таком подхо- 


де возникает возможность применения эффектив- 
ных численных алгоритмов расчета для решения 
каждой задачи в отдельности для всей сборки в це- 
лом и автоматического учета граничных условий. 

С учетом этого основные этапы алгоритма рас- 
чета характеристик ЭГС в дуговом режиме сводят- 
ся к следующему. Пусть в момент времени /=? 0 в 
стационарное решение вносится возмущение по 
одному из параметров математической модели, ко- 
торое приводит к нарушению теплового либо элек- 
трического равновесия. Для расчета переходного 
процесса по исходным ВАХ ТЭП, известной на- 
грузке и температурам электродов определяется 
новое распределение напряжений и токов. 

Далее производится теплофизический расчет 
ЭГС, включающий решение уравнений теплопро- 
водности с учетом распределения тепловыделения 
по высоте, выделения джоулева тепла при прохож- 
дении тока по электродам и коммутационным пе- 
ремычкам, передачи тепла излучением, теплопро- 
водностью цезиевой плазмы и электронным ох- 
лаждением с эмиттера, тепловых потоков с торцов, 
по перемычкам, дистанционаторам и условий от- 
вода тепла от коллектора. 

По результатам расчета температурного поля 
производится перерасчет теплофизических и элек- 
трических характеристик по высоте ЭГС и нахо- 
дится решение для следующего момента времени 
по изложенному выше алгоритму. Процесс расчета 
переходных характеристик заканчивается после 
установления стационарного состояния. 

Анализ существующих методов решения ура- 
внений, представленных в [9], с соответствующими 
граничными условиями показывает, что наиболее 
приемлемым является метод конечных разностей, 
позволяющий свести нелинейную систему диффе- 
ренциальных уравнений в частных производных с 
нелинейными граничными условиями к алгебраи- 
ческой системе уравнений с изменяющимися пра- 
выми частями. Часто при решении подобных задач 
используется явная схема аппроксимации произ- 
водных, обладающая большой наглядностью и про- 
стотой. Однако жесткие ограничения выбора шага 
по времени, присущие явным схемам, в ряде случа- 
ев приводят к неоправданному увеличению време- 
ни счета на ПЭВМ. Поэтому использована эконо- 
мичная схема аппроксимации производных диффе- 
ренциальных уравнений, основанная на методе 
Писмана-Рэкфорда [10], которая позволяет сни- 
зить жесткие требования к выбору шага по времени. 


В соответствии с методом Писмана-Рекфорда 
составлена система линеаризованных алгебраиче- 
ских уравнений для каждого шага по времени и 
каждого узла расчетной схемы і, к вида 


АЛ и 2 - (1 + К КГ + К Г 1 = -К > 


і = 1,..., 2Ѵ(»г(7) + 3) +1, 

- (1 + В 2М К Г 1 + К Т п\ = -< 1/2 ’ (1) 

к = \,...,Ы г -\, п = 0,1,2,... 
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< к =№г(-1 + В у К + В Хк Т" +ГГ< 


С = Л 2 

+ ™. 


-(-і+л,,хг 1/2 + 


и+1/2 
Л ’ 


( 2 ) 


где А и , А гі , А у , В и , В 2І , В ЗІ - коэффициенты конечно- 
разностной аппроксимации уравнения теплопро- 
водности в продольном и радиальном направлениях 
соответственно; /V- число ЭГЭ в ЭГС; т(І) - число 
расчетных точек по длине эмиттера /-ого ЭГЭ; 
/Ѵ г - число точек по радиусу; Т" к - температура в точ- 
ке і, к на временном шаге и; г - шаг по времени. 

Нелинейные функции ю" к , о>"/ 1/2 имеют различ- 
ный вид в узловых точках, расположенных на раз- 
ных элементах конструкции ЭГС (рис. 1). Напри- 
мер, для узлов, расположенных на эмиттерах, в от- 
личие от узлов, расположенных по торцам ЭГЭ, в 
эти функции входят нелинейные члены, описы- 
вающие перенос тепла электронным током и джоу- 
лево тепловыделение. Основным достоинством 
данного метода является то, что стационарное ре- 
шение не зависит от выбора шага по времени г. Для 
решения разностных уравнений теплопроводности 
использован метод переменных направлений, сущ- 
ность которого заключается в следующем. При л=0 
по исходным начальным условиям производится 
расчет правых частей у/" к для всех / и к. Затем мето- 
дом прогонки по коэффициентам А и , А 1Ь Д ѵ вычи- 
сляются температуры на и+1/2 шаге. Полученные 
значения температур используются для расчета 
правых частей (р" к 1/2 . Далее также методом прогонки 
по коэффициентам В ] к , В и , В и рассчитываются 
температуры на л+1 временном шаге. 


Нахождение распределения электрических пара- 
метров ЭГС с учетом нелинейного характера ВАХ 
ТЭП при известных температурах электродов осу- 
ществляется методом Ньютона. На каждом шаге ите- 
рационного процесса ВАХ ТЭП представляются в ви- 
де линейных эквивалентов с напряжениями холосто- 
го хода Е ік и внутренними сопротивлениями К ік [3]. 

Эквивалентная электрическая схема замещения 
участка протяженных электродов приведена на рис. 2. 


ѵ Еіл,к Кщш)Л Ѵ ы ™ Уш + ы 



Рис. 2. Эквивалентная электрическая схема замещения 
участка электродов 


С учетом этого применение конечно-разност- 
ной аппроксимации к уравнениям потенциалов и 
граничным условиям, записанным в [9], позволяет 
представить их в виде 


Так как диагональные элементы в ( 1 ) преобла- 
дают, следовательно, процессы прогонки устойчи- 
вы, а решение разностной задачи существует и 
единственно. К достоинствам метода, достигнутым 
благодаря введению промежуточного слоя, следует 
отнести: второй порядок аппроксимации произ- 
водных, абсолютную устойчивость и расщепление 
исходной задачи на две более простых, решаемых с 
помощью алгоритма скалярной прогонки. 

Описанная процедура при расчете установив- 
шегося температурного поля повторяется либо за- 
данное число раз, либо до достижения требуемой 
точности. 



Рис. 1. Расчетная схема тепловых процессов ЭГС: (х) - расчет- 
ные точки температурного и потенциального полей 


у -у у -у 

г Еі-\,к г Еі,к г Еі,к ѵ Еі+\,к 
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Е(і-1,і),к 
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Е(і,і+1),к 


Ѵ„, -V,, 


' Еі,к ' Сі.к , Е і к (У Еі к Ѵ Сі к ) 


II, 


К, 


-ІЛ= О, 


у _ у у -У у _ у 

ѵ С /-!,* Сі.к ' Сі,к ’ СМ,к + г Еі,к г Сі,к 


в ѵ 


К С (і-1,і),к 2? сц,і+1),к ±ѵ Уі.к 

Е,, к ~ (Ѵ Еіік ~ Ѵ Сі к ) Ѵ СІЛ -Ѵ о 


22 , і- Кумі к 


-/« = 0, 


( 3 ) 


где Ѵ Еік , Ѵ Сік - потенциалы эмиттера и коллектора; 

22с(і-иц - сопротивления электродов между 
узлами (/-1,/) к - ого ЭГЭ; К Гік , К Шк - сопротивле- 
ния утечек тока между электродами и через коллек- 
торную изоляцию соответственно; Ѵ 0 - потенциал 
корпуса; і - номер расчетного сечения к-ото ЭГЭ; 
к - номер ЭГЭ в ЭГС; 

І А = 0 при /'= I , к= 1 и і=2,...,т(к)—\, к=\,...,И, 
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Рис. 3 . Область коммутационной перемычки и ее расчетная схема 


V -V 

т і Ет(М), N 0_ т г тг 

1 А~ 1 К\~ 0 > У Ет(р)+\у — У 0’ 

к п 

где У РЛ - потенциал к - ой перемычки, К РІ к - сопро- 
тивление цилиндрической части перемычки; 
К п - сопротивление токовывода /Ѵ-ого ЭГЭ. 
І в =-І к пщ і= 1, к= 1, І в = 0 при і=1,...,т(к)—1, к= 1,...,ІѴ, 

V -V 

Т ѵ Ст(к-\)+1,к-1 г С\,к . Л , ~ , г 

І в = при г=1, С=2,...,А, 

V. -V 

г Ст (к ),к Ст(к ')+ . /7\7 1 ЛГ1 

І в = — — — при г =т{к),к = 

К ССк 


, ’ ОиСіѴѴѴ *0 . /ДГ Ч 

І в = при і = ш(Л'), 


где Л са - сопротивление утечки тока между кол- 
лекторами А:-ого и А:+ 1 -ого ЭГЭ. 

Для коммутационных перемычек (рис. 3) мож- 
но записать 

^ Ет(к\к — Ѵр.к ^т(к)+1,к ~ ( ^Е\,к+\ ~ Г Р,к ) 


Г -Г Г -V 

_І_ ' ги+1 ' Р,к _ г Р,к г Ст(к)+\,к _ д 


V -V V -V 

ѵ Ст(к),к г Ст(к)+1,к г Р,к г Ст(к)+1,к 


V -V 

г Ст(к)+\,к ’ С1.І-+1 


^ Ст(к)+1,к 

при і = т(к) + 1,к = 1,...,іѴ-1, 


где К т - сопротивление дисковой перемычки. 

Распределения токов, потенциалов и напряже- 
ний с учетом токов утечек между электродами, эл- 
ементами и на корпус находятся путем решения 


линеаризованных уравнений (3, 4) с заданными 
граничными условиями. В случае отсутствия соот- 
ветствия между электрическими характеристика- 
ми, полученными в результате решения системы 
уравнений и характеристиками, определенными по 
ВАХ ТЭП, вычисляются новые значения электри- 
ческих параметров линейных эквивалентов харак- 
теристик термоэмиссионных преобразователей. 

Итерационный процесс решения системы уравне- 
ний (3, 4) с внесенными поправками по Е ік и Д к про- 
должается до установления требуемого соответствия. 

Основным достоинством описанного метода 
расчета электрических характеристик ЭГС являет- 
ся то, что при любом изменении параметров элек- 
трической схемы автоматически пересчитывается 
распределение потенциалов и токов во всех ЭГЭ. 

Таким образом, разработанные методы расчета 
характеристик многоэлементных ЭГС, позволяют: 

• определять распределение температур в вакуум- 
ном и диффузионном режимах работы незави- 
симо от величины генерируемого тока. В этих 
случаях при расчетах ВАХ ЭГС температурное 
поле считается неизменным; 

• решать задачу по определению согласованного 
температурного и электрического полей в раз- 
рядном режиме работы с учетом электронного 
охлаждения электродов, джоулева тепловыделе- 
ния при протекании тока по электродам и ком- 
мутационным перемычкам; 

• получать серии ожидаемых ВАХ для всех режи- 
мов испытаний в широком диапазоне измене- 
ния давления паров цезия, уровней удельного 
тепловыделения, расхода и входной температу- 
ры теплоносителя; 

• моделировать путем изменения параметров 
расчетной схемы возникновение и развитие 
аварийных ситуаций в ЭГС. 
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